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Abstract. We study the metal- insulator transition on a three dimensional quantum per- 
colation model by analyzing energy level statistics. The quantum percolation threshold pq, 
which is larger than the classical percolation threshold pc, becomes smaller when the time 
reversal symmetry (TRS) is broken, i.e. Pq(with TRS) > pq (without TRS) > Pc- It is shown 
that critical exponents are consistent with the result of the Anderson transition, suggesting 
that the quantum percolation problem can be classified into the same universality classes of 
the Anderson transition. The shape of level statistics at the critical point is also reported. 
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1 Introduction 

Due to its simple but non-trivial nature, percolation problems have been attracting a 
lot of attention . Especially interesting is its application to the transport properties 
in three-dinemsional disordered solids. However, at low temperature, the quantum 
interference effects should be seriously taken into account and the classical percolation 
picture is insufficient. Thus the study of the quantum percolation problem becomes a 
very important subject. 

It is well known that disordered systems show the Anderson metal- insulator tran- 
sition. This is classified into three universality classes (orthogonal, unitary, and sym- 
plectic universality classes) by time-reversal and spin-rotation symmetry. This con- 
jecture is well confirmed for the Anderson tight binding model [^. Then it is natural 
to ask whether the quantum percolation problem can be also classified into the same 
universality classes of the Anderson transition. If yes, the critical exponent for the 
divergence of the localization length in the quantum percolation model is the same as 
in the Anderson model, namely, v — 1.57 ± 0.02 in the presence of the time reversal 
symmetry (TRS) ||, ||, and = 1.43 ± 0.02 in the absence of it g. However, the 
estimates of v in the quantum percolation problem reported so far are at variance. 
For example, renormalization group analyses of the quantum percolation problem 
give u = 2.1 |, |, 1/ = 1.9 ± 0.5 §, and v = 1.86 ± 0.02 §, weU above the value 
in the Anderson model. From the Thouless number analysis, it is estimated to be 
— 1.95 ± 0.12^, |j, again larger than in the Anderson model. On the other hand, 
the transfer matrix method for the network model gives v = 0.75 ± 0.1 llQ], and the 
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analytic estimate of the transmission coefficient gives v = 0.38±0.07[|ll), significantly 
smaller than in the Anderson model. Recently, using the finite size scaling behavior 
of level statistics Berkovits and Avishai estimated to be 1.35 ± 0.1 not far 
from the value of the Anderson model. Inspired by their work, we perform extensive 
numerical study of the level statistics on the percolation cluster. We also study the 
effect of breaking of TRS. 



2 Model and method 

To describe the three dimensional (3D) quantum bond percolation model, we consider 
the following tight-binding Hamiltonian 

H = Y,{U,a\a,^^A.c) (1) 

where (if) denotes nearest neighbors. The transfer energy is defined as 

^ f exp(i0jj) (for connected bond) 

[ (for disconnected bond) 

Bonds are randomly connected with probabilities p. Oij is the Peierls phase due to 
magnetic field. When all the Peierls phases are set to 0, the Hamiltonian is time 
reversal symmetric, and we call it TRS model. When the phases are not vanishing, 
the Hamiltonian is generally not time reversal symmetric. We set — tt < dij < tt 
randomly, and call this situation non-TRS model hereafter. The underlying lattice 
is a three-dimensional cube of length L with periodic boundary conditions. For each 
realization of bond structure, we first find maximally percolating cluster and then 
we diagonalize it by Lanczos method. In order to gain good statistics, more than 10^ 
eigenvalues are calculated, so the number of realizations of random bond configuration 
are TV = 580, 300, 175 and 110 for sample sizes 12^, 15^, 18^ and 21^. We take 

levels in the region 0.2 < \E\ < 0.8 where density of states is rather smooth [ p^ . 
Then the distribution function P{s) of adjacent level spacings s is calculated. In the 
limit of large system size L — > oo, the level spacing distribution P{s) is expected to 
be described by Poisson distribution P(s) = exp (— s) in the insulator region. In the 
metallic regime, P{s) is well described by the Wigner surmise, P(s) oc exp (— ^^s^) 
(/3 = 1 in the presence of TRS and /? = 2, otherwise). 

In order to obtain the critical value of the probability pq and critical exponent v, 
we define I{s) and A{p, L) as 

Jo Iwisisjds-J^ Ip{s)ds Jo 

which characterizes the transition from Poisson to Wigner. sq is set to 1.2. Noting 
that ^(p) diverges as 
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Fig. 1 Level spacing distribution for TRS (a) and non-TRS (b) models at the critical point. 
The system size is L = 21 for example, but the form is almost independet of L. 

near the critical probability pq, A(p, L) is expected to behave as 

A(p, L) = f[L/ap)] = ao + a,[Lip - p^Y]'/-" + a^[L{p - p^ff/'' + • • • . (5) 
Fitting the data to this expression, we estimate pq and v. 

3 Result 

The results are summarized in Table I. When magnetic fields are applied, the critical 
point Pq of non-TRS model becomes smaller than that in the presence TRS, suggesting 
that this transition has a nature of the Anderson transition. Note that Pq is much 
larger than the classical bond percolation threshold Pc ~ 0.249. This is why taking all 
energy levels including small clusters fl^ and taking only those in the lagest cluster 
make small difference |l^. Both TRS and non-TRS models show v consistent with 
the value of the Anderson transition . 

We then compare the shape of level statistics at the critical point pq(Fig. 1). It is 
clearly seen that the critical level statistics is sensitive to the breaking of TRS. 

In the large s region, the level statistics is described by the sub-Poisson form, 
(X exp (— Ks) (Fig. 2). In the Anderson model, k, — 1.9 ± 0.1 for the orthogonal 
symmetry [|l^, and 1.87 for unitary symmetry |l^. In the present quantum percolation 
model, K = 1.77±0.10 in the presence of TRS, and for non-TRS model n = 1.87±0.03. 

The overall form of the critical level spacing distribution deviates from those in 
the Anderson model with periodic boundary condition (b.c.) and the fixed b.c. [|l5j. 
Detailed will be reported elsewhere. 
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Table 1 The critical point pq and the critical exponent v with their standard deviations. 
The values of the TRS model and the non-TRS model are only for percolating cluster, while 
TRS (all) is obtained from all energy levels. TRS (all,BA) is the result reported in ref. [0|. 



system 


size 








V 




TRS (all,BA) 


L = 


7,9,11,13,15 


0.33 ± .01 


1.35 ± 


.10 


TRS (all) 


L = 


12,15,18,21 


0.321 ± 


.001 


1.45 ± 


.07 


TRS 


L = 


12,15,18,21 


0.324 ± 


.001 


1.42 ± 


.07 


TRS 


L = 


15,18,21 


0.324 ± 


.001 


1.45 ± 


.11 


non-TRS 


L = 


12,15,18,21 


0.309 ± 


.001 


1.13± 


.05 


non-TRS 


L = 


15,18,21 


0.308 ± 


.001 


1.25 ± 


.08 
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